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Abstract 

An analysis of magnetic anisotropy in RC05 compounds is performed with 
account of screening of ion point charges by conduction electrons. A crucial 
role of non-uniform distribution of screening electrons (the terms containing 
derivatives of charge density) is demonstrated. Influence of anisotropy of 
screening, that is connected with the anisotropy of the Fermi surface, on the 
magnetic anisotropy sign is discussed. 
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It is believed now that the main mechanism of occurrence of strong magnetic anisotropy 
(MA) in the rare-earth (RE) based intermetallic systems is the crystal field one The 
simplest point-charge model leads frequently to difficulties and contradictions with experi- 
mental data. For example, the calculated MA constant K\ in RC05 compounds turns out 
to be very large and have an incorrect sign ||. Similar difficulties occur for the R 2 Fei4B 
systems 0. Thus screening of the crystal field should play an important role. This screen- 
ing can be taken into account by introducing the effective ion charge Q* which can differ 
considerably from the bare ion charge. Shielding of crystal fields in ionic solids || turns out 
to be insufficient even for qualitative explaining experimental data. 

Another approach to the MA problem is based on the first-principle band-structure 
calculations. Generally, RE systems are a difficult case for standard band theory methods, 
so that one has to use some approximations and model representations (e.g., atomic sphere 
approximation in early papers, inclusion of /-states in the core states). Modern calculations 
IIU treat mainly the anisotropic contribution of the conduction electrons in the atomic 
sphere with the centre at the RE /-ion, and the rest "lattice" contribution is obtained in 
the point-charge approximation with screened charge densities of other cells, which do not 
overlap with the /-shell. It should be noted that the atomic and lattice contribution are in 
fact not independent, and their partial cancellation takes place. The local- spin-density full- 
potential approaches, which are based on total energy calculations |J, can provide correct 
orders of magnitude of anisotropy energy, but quantitative agreement with experimental 
data is not satisfactory. 

Thus simple qualitative considerations which start from the physically transparent point- 
charge picture, but introduce effective distance-dependent ion charges, seem still be useful. 
In the present paper we analyze the first MA constant with account of charge screening 
by conduction electrons in a metal. As demonstrated in Ref. effective charges depend 
strongly on the concrete form of screening electron density. In particular, the Friedel oscil- 
lations of charge density can play an important role. Here we discuss the screening in more 
details, in particular with account of its anisotropic character. 
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The MA constants are determined from angle dependence of the energy of magnetic ions 
in the crystal field 

5S cf = -K lC os 2 6 + ... (1) 

We consider the magnetic ion at the point r = in the crystal field of the surrounding 
charges. The contribution to the crystal-field potential from the ion with the bare charge 
Qq at the point R can be represented as 

= 0o + q, ( R- r) 

|R — r| 

where Q e i(TV) is the conduction-electron screening charge. After expanding @ in r up to 
r 2 , the expression for K\ can be written in the form (cf. Refs. 

K{ c = 3e 2 A(r 2 f )ajJ(J - 1/2) (3) 

Here (r 2 ) is the average square of the /-shell radius, J is the total angular momentum of 
the RE ion, aj is the Stevens factor, 

A^ W^-' , (4) 

R n 

the sum goes over the lattice, #r, is the polar angle of the vector R, Q*(R) are the corre- 
sponding screened ion charges. Note the difference of our consideration from Ref. ||, where 
the screening charge in the sphere with the centre at the magnetic RE ion was calculated. 
In the spherical charge density approximation we have [|7| 

Q*(R) = Q + Q e i(R) - ^nR 3 [Z(R) - RZ'(R)] (5) 

where Q e i(R) is the conduction electron charge inside the sphere with the centre at the point 
charge and radius R, 

Q el (R)=A7r J% 2 dpZ(p) (6) 

Z(R) is the charge density, Q' el (R) = 4irr 2 Z(R), the system of units where the electron 
charge e = — 1 being used. We see that Q*(R) includes explicitly, besides the total charge 



Qei(R), also the charge density Z(R) and its derivative Z'(R). Such terms do not occur in 
the calculations H where the "lattice" charge density does not overlap with the /-shell. 
Note that higher-order anisotropy constants are obtained after keeping next-order terms in 
r and include higher-order derivatives of Z(R). 

To obtain the value of Q*(H), one has to investigate the charge screening for a concrete 
electronic spectrum. In Ref. 0, the one-centre screening problem was considered within a 
simple model of free conduction electrons (E = k 2 /2) in the rectangular potential well which 
is induced by single impurity and has the width d and depth Eq = k\j2 |§. This model 
enables one to calculate the charge distribution of screening conduction electrons in terms 
of the spherical Bessel and Neumann functions ji(kr) and ni(kr) (r > d) and the scattering 
phase shifts rji- The value of k is determined for given k F and d from the Friedel sum rule 

c\ OO 

Q = -J2(2l+l) m (k F ) (7) 
71 1=0 

The parameter d should be determined by the geometry of the lattice near the impurity. 
In Ref. ||, where impurities in the Ag host were considered, d was chosen to be equal to 
the Wigner-Seitz radius, so that k F d = 2. The results of the calculations for kpd = 2 and 
k F d = 3 are presented in Figs. 1,2. At Qq = 1 Eq.(0) yields kod = 1.46 and hod = 1.235 
respectively. 

The model discussed is more adequate for impurities which induce strong disturbance 
of charge density (e.g., in hydrogen-containing RE systems fHj). Of course, the choice of d 
may be different and more complicated models are required for a regular lattice of screened 
charges where interference of screening charge clouds from different centres takes place. In 
metals, the value of ion charge Qq can be put to the charge in some sphere (e.g., charge 
transfers in the atomic sphere were considered in Ref. H) and does not necessarily coincide 
with the nominal free ion value. Besides that, the dependence Q*(R) in the lattice can 
become anisotropic. Nevertheless, we use Figs. 1-2 for a qualitative discussion. 

One can see from Fig. 2 that at R < d, except for the case of very small R where Q*{R) 
slightly decreases, the derivative term in (f|) results in that Q*{R) grows (despite an increase 
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of \Q e i(R)\). For R ~ d, where Z' is maximum, the non-uniform distribution of electron 
density leads to that the effective ion charge Q* is positive and exceeds considerably its bare 
value Qo (Qo = 1 in our case). At the same time, with further increasing R the situation 
changes drastically: Z' decreases and becomes negative, so that "overscreening" of the ion 
charge occurs. At large distances Q* tends to zero, but considerable oscillations of the 
effective charge sign take place, which attenuate rather slowly. It should be stressed that 
the oscillations are due mainly to the derivative term. On the other hand, the quantity 
Qo + Qei(R) monotonically decreases up to R ~ d and then tends to zero very rapidly, 
oscillations being very weak (Fig.l). This agrees with the fact that band calculations yield 
usually small values of charge transfers in atomic spheres (see, e.g., 0). The oscillation 
period and position of maxima and minima of Q*{R) turn out to be weakly sensitive to the 
values of Qo an d d, but are determined mainly by the parameter kp. 

Now we analyze concrete geometry of the RC05 lattice (CaCus structure, Fig. 3) with 
c ~ 4A, a ~ 5A. The Co ions have two positions in two different types of hexagonal layers: 
2c (Col) sites in the layers containing R atoms with R = a/\f?> ~ 0.57a and 3g (Coll) sites 
in layers with no R atoms with R = |v 'a 2 + c 2 ~ 0.64a. Although positions of all R ions 
are equivalent, the contributions to the crystal field at a given R site from the charges in 
the same plane (sites RI, R = a) and in neighbor planes (sites RII, R = c) can be different, 
since the effective charge is a function of the distance between two sites. Then we can write 
down (cf. [|]) 

A = 6a- 3 (l6 (1 2 ^~^ /2 gU - ^ ,2 Qhoi + \y-"Qlu - Qri) (8) 

where y = c/a ~ 0.8. The contributions of in-plane (I) and out-of-plane (II) neighbors of R 
site have different signs. However, unlike the pure rare earth-metals with the hep structure 
]]]], the small geometrical factor ^8/3 — c/a = 1.633 — c/a ~ 0.05 does not occur, so that 
the calculated value of K\ turns out two orders of magnitude larger. 

The experimental data on the c/a ration and first magnetic anisotropy constant of the 
RCo5 +a: compounds at low temperatures are presented in Table 1. A monotonic decrease of 
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c/a in the RE series takes place (with exception of tetravalent Ce). 

Table 1. The total angular momenta J, Stevens factors aj, average squares of the f- 
shell radius (r"j) (atomic units) for free R ions; c/a ratios and the experimental values of 
K\ (K/RE ion) in RCos +x systems [for light R according to fI^J12|| and for heavy R (the 
contribution of RE sublattice) according to Ref. jl3"| ]. The corresponding values of Q* Co are 
calculated by using (|]), (||). 



RCo5 +x 


CeCo 5 PrCo 5 NdCo 5 SmCo 5 TbCo 5 .i 


DyCo 5 . 2 


H0C05.2 


ErCo 6 


J 


5/2 4 9/2 5/2 6 


15/2 


8 


15/2 


aj x 100 


-5.7 -3.4 -7.1 4.1 -1.0 


-0.63 


-0.22 


0.25 


(rj) 


1.20 1.09 1.00 0.88 0.76 


0.73 


0.69 


0.67 


c/a 0.817 0.797 0.796 0.795 0.803 0.810 0.817 0.821 




-61 -44 -220 190 -96 


-211 


-203 


80 


Qco 


-0.007 -0.003 -0.006 -0.04 -0.02 


-0.03 


-0.09 


-0.04 



One can see that the experimental values for heavy rare earths and Sm are obtained 
for Aa 3 ~ 1. Taking into account only the contributions of Co ions in (|8]) and assuming 
Qhoi = Qcoii = Qco we obtain 



A 



23.4^* 



(9) 



Thus, for Q* Co ~ 1, K\ has incorrect sign and is very large in absolute value (of order of 
3000K). It should be noted that such values of K\ are in fact not self-consistent: they would 
destroy the Russel- Saunders coupling and quench total momenta. In light rare earths (Ce, 
Pr, Nd and Sm) the momenta, as obtained from both neutron scattering and magnetization 
measurements |L4j , are indeed considerably suppressed. Ce ions are supposed to be tetrava- 
lent; for other light rare earths, effects of the strong crystal fields and exchange interactions 
can play a role. On the other hand, in Tb, Dy, Ho and Er the saturation momenta are close 



to their free ion values 15 



To compensate large numerical factor in (^), one has to put Q* Co to be very small in 
absolute value and negative (Table 1). One cannot rule out that the distance between Co 
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ions corresponds to negative values of Q* with \Q*\ < 1 in Fig.2, but such a situation 
is rather unusual. Thus the simple model with equal Co charges can hardly explain the 
observed sign and value of K±, and a more detailed treatment of screening is needed. 

An assumption that main role in the crystal field formation belongs to R ions and Qq q 
can be assumed zero was made in Refs. |T5|,|I2|] . The effective charges for R ions, Q^, could 
be supposed to equal about 1, as well as in pure rare earths where this parameter varies 
between 1.1 and 1.4 0. Then we would obtain the correct sign and order of magnitude for 
K\ since the contribution of RII ions dominates. However, the distances between RE ions 
in RC05, c and a, are considerably larger than in pure rare earths (about 3.5A), so that the 
values Q\ can be much smaller. Besides that, the variation of effective charge in the RC05 
series is much stronger than for pure rare earths. Such a variation can be related to the 
contribution of Coll sites since the factor at the first term in the brackets of (§) depends 
appreciably on y = c/a. 

It is instructive to take into account the anisotropy of screening of Co charges that is con- 
nected with the anisotropy of the Fermi surface in the hexagonal lattice. Then the screening 
is determined by the wavevector of the Fermi surface in the corresponding direction. The 
effective values of fc_p in the hexagonal planes can be supposed smaller than in the direction 
to Coll sites. Then the charge of the Coll ions is screened with distance slower than for Col 
ions and can dominate in Ki, despite the larger distance and numerical factor in (H). This 
circumstance may lead also to additional (besides pure geometrical) strong dependence of 
the anisotropy on y. It should be noted that the calculation for SmCos || yields for charge 
transfers in the Co atomic spheres (as well as for Col and Coll magnetic moments) different 
values: g C oi - 0, g C oii - -0.03. 

Experimental values of K 1 from Table 1 can be obtained, e.g., for Qq oU — 4Q^ oI ~ 0.3, 

= 0, or for Qq oU — 3Qc i — 0-2, ~ 1 (in the latter case, the sign of K\ is determined 
by R contribution which is partially cancelled by the Co contribution). Thus occurrence of 
very small values of Q* Co by accidental reasons is not required in such a consideration. 

To conclude, an account of distance-dependent screening of ion charges by conduction 
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electrons can provide a reasonable explanation of the first MA constant value in RC05 
compounds, the non-uniform distribution of screening electrons (the derivative term in the 
expression for the effective charge (|5])) being of crucial importance. Our approach does not 
introduce artificial separation into atomic spheres. On the other hand, it has a number 
of other drawbacks. In particular, it treats independent charge centres and is not self- 
consistent: perturbation of charge density by R ion itself (see, e.g., Ref. |§) is not taken into 
account. Thus a synthesis of our consideration with approaches, that are based on the real 
electronic structure calculations, would be of interest for further insight into the problem of 
magnetic anisotropy in the rare-earth systems. 
Figure captions 

Fig.l. The distance dependence of the sum of the bare ion charge and electron screening 
charge, Q + Q ei (R), for k F d = 2 (solid line) and kpd = 3 (dashed line) with Q = 1. 

Fig. 2. The distance dependence of the effective charge Q*{R) © for the same parameter 
values as in Fig.l. 

Fig. 3. The local environment of the R site in the crystal structure of RC05 compounds. 
Large open circles denote R sites, small open circles Col (2c) sites, small shaded circles Coll 
(3g) sites. 
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